Spatio-temporal spectral analysis of a forced cylinder wake 
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The wake of a circular cylinder performing rotary oscillations is studied using hydrodynamic 
tunnel experiments at Re — 100. Two-dimensional particle image velocimetry on the mid-plane 
perpendicular to the axis of cylinder is used to characterize the spatial development of the flow and 
its stability properties. The lock-in phenomenon that determines the boundaries between regions of 
the forcing parameter space were the wake is globally unstable or convectively unstable (see [1 for a 
review) is scrutinized using the experimental data. A novel method based on the analysis of power 
density spectra of the flow allows us to give a detailed description of the forced wake, shedding light 
on the energy distribution in the different frequency components and in particular on a cascade-like 
mechanism evidenced for a high amplitude of the forcing oscillation. In addition, a calculation of 
the drag from the velocity field is performed, allowing us to relate the resulting force on the body 
to the wake properties. 

PACS numbers: 47.15.Tr 47.80.Cb 47.20.Ft 



I. INTRODUCTION 



The simple geometry and the complex behavior of the 
flow around a circular cylinder at low Reynolds numbers 
{Re = DUq/v < 180, where D is the diameter of the 
cylinder, Uq the free-stream velocity and v the kinematic 
viscosity of the fluid), makes it a prototypical 2D wake 
flow. The well-known Benard-von Karman (BvK) vortex 
street [2] |3] results from the destabilization of the steady 
flow in the wake of the cylinder, driven by the periodic 
shedding of opposite-signed vortices, that occurs above 
the threshold Rec ~ 47 (see e.g. [4, 5j). Often used as 
a model for shear flow instabilities (see e.g. [6HTQ]), the 
cylinder wake gives a framework to study the distinctive 
features of spatially developing flows. In particular, for 
Re > Rec the velocity field in the whole flow domain os- 
cillates with the same global frequency and its harmonics 
and, because the oscillation is spatially evolving, it can 
be characterized through the evolution of its envelope or 
global mode of the instability as a function of the flow 
parameters [TT| [T2] . 

The cylinder wake has been also widely used to test meth- 
ods of flow control using dynamic actuation, either in an 
open- loop sense [13j or with a feedback loop closed by the 
signal taken by a flow sensor [14J and within an optimal 
control scheme [15j. A large subset of the existing liter- 
ature about control strategies on the cylinder wake con- 
cerns the use of imposed oscillations, the case of in-line 
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oscillations [16j being intimately related to the problem 
of vortex-induced vibration (see extensive review by [17]). 
For rotational oscillations, from the first visualizations of 
[T8] and the experimental work of [19j, to the numerical 
[2QH22] and experimental works [23j that have followed, 
it has been shown in particular that the imposed rota- 
tional oscillation can significantly modify the geometry 
of the cylinder wake and hence the drag coefficient. The 
stability properties of the forced wake have been studied 
by [24J using experimental measurements of the velocity 
field. They described in terms of the (/, A) parameter 
space (the frequency / and the amplitude A of the oscil- 
lations) the two qualitatively different states that arise 
in the wake: a spatial mode dominated by the BvK vor- 
tex street, meaning that the wake is globally unstable; 
and the so-called lock-in regime, where the frequency is 
imposed by the forcing in the near wake, the amplitude 
of the oscillation rapidly decaying downstream with the 
characteristics of a convective instability. 
The first goal of the present work is to refine the de- 
scription of the transitions between the locked regime 
and the global instability, which have been described in 
previous works on this setup (see [Ij for a review), espe- 
cially for forcing frequencies lower than the natural fre- 
quency. Particle image velocimetry (PIV) measurements 
allow us to study the spatial development of the forced 
flow, in particular the modification of the global mode 
represented by the velocity fluctuation envelope. Scaling 
laws in the forcing parameter space that universally de- 
scribe the evolution of the global modes in the vicinity of 
the critical lines separating the lock-in and the globally 
unstable regions are derived from this data. A spectral 
analysis of the velocity fluctuations gives an alternative 
procedure to define the critical lines and confirms previ- 
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FIG. 1: (Color) Top: Experimental setup and definition of 
the forcing parameters. Bottom: Vorticity field calculated 
from PIV measurements for the non-forced case showing the 
Benard-von Karman vortex street. 

ous linear stability studies. The drag forces for each case 
are estimated from the PIV measurements using the so- 
called flux equation (see [25j). These results allow us to 
relate the resulting force on the body to its wake proper- 
ties. 
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FIG. 2: (Color) Forcing parameter space. The points repre- 
sent experiments where: ■ no lock-in is observed; ^ lock-in 
is observed. The blue line is the lock-in threshold estimated 
from the present experiments. The red line stands for the ob- 
servations of [24] . Striped region indicates a maximum drag 
reported by [26] . 



II. EXPERIMENTAL SETUP AND 
PARAMETERS 

The experimental setup shown in Fig. [T] is the same 
used by [23j. A circular cylinder of diameter D = 5mm 
and span of ^ 20D is placed in a low-speed hydrody- 
namic tunnel with a 100 x 100mm cross-section. The 
cylinder span thus practically covers the whole height of 
the tunnel and the ratio of the cylinder diameter to the 
tunnel section width is of 1/20. The cylinder can perform 
rotational oscillations driven by a stepper-motor placed 
below the platform of the test section on a submerged 
'technical section' of the tunnel, imposing a controlled 
forcing on the flow. The Reynolds number is set to 100, 
defining it using the up-stream velocity in the center of 
the tunnel as Uq. The measured natural vortex shedding 
frequency was /o = 0.63 Hz, so the Strouhal number 
St = foD/Uo is approximately 0.15. Quantitative mea- 
surements were performed using 2D particle image ve- 
locimetry (PIV) on a horizontal plane placed at mid-span 
of the cylinder (see the vorticity field for the non-forced 
case which corresponds to the well known Benard-von 
Karman (BvK) vortex street in Fig. [T]). Image acquisi- 
tion and PIV calculation were done using a LaVision® 
system composed of an ImagerPro 1600 x 1200 CCD cam- 
era with a 12-bit dynamic range recording double-frame 
images at llHz and a two rod Nd:YAG (15m J) pulsed 
laser synchronized by a customized PC using LaVision 
DaVis 7.1 software. Laser sheet width was about 1 mm 
in the whole 100 x 80 mm imaging region. The time lapse 
between the two frames of each image pair used for PIV 
was set to At = 12ms and sets of 500 snapshots give a 



frequency resolution of A/ = 0.02Hz. 

The rotational oscillation of the cylinder is prescribed 
by a forcing function of frequency / and amplitude 
that can be written as 0{t) = cos(27r/t), which allows 
the forcing to be unequivocally described using two inde- 
pendent non-dimensional parameters as did by [TS]: the 
forcing amplitude A = uoma^/Uo, where uoma^ = DirfOo 
is the maximal azimuthal velocity of the rotational oscil- 
lation; and the ratio f/fo. We explore 100 forcing cases 
and they are represented in Fig. |2j We define Ac as the 
value of the amplitude for the transition between lock-in 
and non-locked state (critical lines). These values are ob- 
tained from spectral analysis developed on section [V| and 
Fig [2] compares them with previous works based on linear 
stability from [24j at Re = 150. In order to compare the 
lock-in region with drag force estimation (section IV), 
we include in the figure the maximum drag found in^26j 
numerical work for Re = 200 



III. REVISITING VORTEX PATTERNS AND 
GLOBAL MODES IN THE WAKE 

A first step to characterize a regime that arises un- 
der forcing is by identifying its vortex patterns, as ini- 
tially done in the early works of [16j. In Fig. [sj we show 
the wake patterns for two different forcing frequencies 
and various amplitudes. In the left column of Fig. [3) 
// = 0.69/o and Ac ~ 0.3. We observe at low amplitude 
(a) A = 0.10 that the flow is not entirely synchronized 
with the natural frequency but the formation region is 
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// = 0.69/0 



// = 2.40/0 




FIG. 3: (Color) Vortex structures from vorticity contours for 
ff = 0.69/o (Left) and // = 2.40/o (Right): (a) A = 0.10 < 
Ac, {h)A = 0.40 > Ac, (c) A = 1.00 > A^ (d) A = 3.2, (e) 
A = 0.40 < Ac, (f) A = 0.80 < Ac, (g) A = 1.60 > Ac, (h) 
A = 4.00. [See text] 



modulated by the forcing frequency. Increasing the am- 
plitude to (b) A = 0.40, the flow is locked on the forcing 
frequency and vortices are shed in a shorter distance. 
The vortex structure remains the same as the non-forced 
case as each half a cycle a vortex is fed into the wake. 
For A = 1.00 the convective character of the instability 
is clearly evidenced by the intense vortices shed in the 
vicinity. For A > 1 two vortex rows form as the vortex 
cores move away from the cylinder. Once this pattern is 
attained, under further augmentation of the amplitude 
forcing, the vortex structures start splitting in a cascade- 
like pattern as can be seen in Fig. 3] (d). The spectral 
analysis performed below gives more insight on this phe- 
nomenon. The right column of Fig. [3] corresponds to the 
case of ff = 2.40/o and Ac ~ 1, and the evolution of 
the wake pattern under forcing develops in a different 
way. For A = 0.40 (Fig. [sj^e)), the wake aspect does not 
differ from natural BvK vortex shedding except that in 
the near wake the forcing wavelength slightly modulates 
the vortex formation. As the forcing amplitude increases 
the shedding distance from the cylinder also grows: in 
Fig. [3|f) it is - lOD, for A = 0.80 < A^ i.e. approx- 
imately twice of what is observed ai A = 0.40. Indeed, 
vortices of the same sign coalesce in the near wake con- 
stituting a long formation region until they are shed in a 
BvK-like pattern farther downstream. Once the lock-in 
threshold amplitude (Ac ~ 1) is exceeded, the coales- 
cence pattern prevails and two vortex sheets are formed 
stabilizing the wake (see Fig.jsj^g)). Higher forcing ampli- 
tudes, for instance A = 4 in Fig. |3jh), destabilize again 
the wake and a BvK-like pattern reappears in the far 
wake. 



A. Global mode shape 

Wake flows can be analyzed as a propagating wave 
with an amplitude (determined by the fluctuating com- 
ponent of velocity) that grows from the origin, reaches 
a maximum and decays afterwards. The spatial enve- 
lope of this coherent oscillation gives the amplitude of 
the so-called global mode, for which the dominant con- 
tribution is given by the first harmonics. Previous works 
[11, 12j studied scaling laws for the global mode in wake 
fiows near the threshold Rcc- A typical contour is pre- 
sented on Fig. EFa) where its maximum amplitude, a^ax 
at (xmax 7 ^ma Jy coordinates, is highlighted as it repre- 
sents an important parameter for scaling . A synthesis for 
the global mode properties is represented by its envelope. 
Fig. |4ja, inset), that corresponds to the position of x^ax 
the maximum amplitude of the mode at i/ — Umax- 

The 

authors proposed that both amax and Xmax follow scaling 
laws such as amax ^ £^ and Xmax ^ £^ in the vicinity of 
the bifurcation. As e represents a control parameter that 
measures the distance to the threshold they determined 
that amax - {Re - Re^'^l^ and Xmax - {Re - i^ec)"^/^. 

The existence of global modes is a consequence of the 
existence of an absolute instability region, which is mod- 
ified by the forcing. The region disappears after lock-in 
is attained. For f > fo this is discussed in [24J and re- 
visited in the present experiment. Fig. |4jc) shows the 
vanishing of the global mode as its envelope broadens 
and its maximum diminishes. Once the lock-in threshold 
amplitude is exceeded, the mean energy of the fiuctu- 
ations Urms — Jo ^x^t/T is driven exclusively by the 
forcing, the curve strongly decreases from the cylinder to 
the wake. This observed envelope is expected to decay 
exponentially with the form Ae^^^, where ^ has the same 
critical behavior as Xmax but with |e|, as it was studied 
in early works by [27j . 

The picture is less clear for the cases where ff < /o, 
shown in Fig.|4jd). While small forcing amplitudes deter- 
mine the global mode maximum to decrease, for A > 0.15 
the envelope increases, and it does it significantly after 
the lock-in threshold. 



B. Scaling laws 

The evolutions of the maximum of the global mode 
envelope amax and its position Xmax can be used to de- 
fine the critical lines that bound the lock-in state in the 
(//,A) space. For a given frequency, the critical am- 
plitude A = Ac can be determined by inspecting the 
evolution of the ratio amax/xmax- in Fig. |4jb, inset). 
As ff is fixed, the effective control parameter changes 
with A, so we expect that amax ^ {A — A^)^!^ and 
^max {A - AcY^I^ which means that amaxl^max 
should behave linearly with {A — Ac). We observe that 
the scaling holds near the critical value, while far from 
the critical lines it is modified by higher non-linearities. 
We scale amax/ ^max and the forcing amplitude A to the 
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distance to threshold with the critical value for each forc- 
ing frequency {A — Ac)/ Ac. Fig. ^h) resumes the three 
cases and a single line represents the linear behavior as 
dmax/xmax ^ for the coutrol parameter {A — Ac) /Ac 
as it attains the threshold. 



IV. DRAG ESTIMATION 
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Previous works on this subject have addressed the 
problem of estimating forces from velocity fields only, em- 
ploying different methods to include the contribution of 
the pressure field in the momentum balance in a control 
volume equation (a framework originally developed by 
[28J): 
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udv 



{-pi + T) • fids 



(1) 



where y is a control volume, S its boundary, p is the 
fluid density, p is the pressure field, I the unit tensor 
and T = ii{Vu + V^u) is the viscous stress tensor. 
The pressure field can be obtained either by means of 
the Poisson equation (see e.g. [29j), or integrating the 
Navier-Stokes (NS) equation along the control surface 
[30 1 . Considering for the pressure p along a 5— curve, 
vi^) = P{s ~ ds) + Vp • ds, the latter idea was further 
refined by [31 J who proposed using the NS equation only 
in the wake region, while adopting the Bernoulli equation 
in the surrounding slowly-evolving potential flow region. 
Thus, p{s) = ^ -\- Po — ^ p\u\'^ which reduces the numer- 
ical error introduced by derivations. 

Another approach has also been used to evaluate the 
force using only velocity fields and their derivatives (see 
e.g. [321 [33] ). It makes use of the the identities (see also 
[34]): 
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where <l> = [{xu)I — xu]^ so that Eq. ([T]) leads to ex- 
pressions where p does not appear explicitly. For a 2D 
problem, the mean flow drag forces becomes: 



FIG. 4: (Color) (a) Global mode spatial structure (from [11'): 
isocontours oi < Ux > for the non-forced flow around a cylin- 
der and shape of the global mode envelope evaluated at the 
maximum (blue region) y — i/max in the inset (the solid line 
is from the present experiments and the points form [TTIV 
From the present experiments: envelope of the global mode 
given by the amplitude of the Ux fluctuations, as a function 
of the downstream distance x to the cylinder and for different 
values of the forcing amplitudes A for // = 2.40/o (c) and 
// = 0.69/o (d); and (b) scaling laws for the spatial envelopes 
for // > fo. Estimation of Ac from amax/xmax versus A in 
the inset and scaled result for amax/xmax in (b). 



(F) = (f ft - {- {u'^)l- (uu) - {u{x X uj)) (4) 
Js 2 

+ [(f • V • (T))I - fV • (T)] + (T) }dS 

where brackets indicate a time averaging procedure. We 
tested both the mixed Bernoulli-NS scheme, and the 
expression from Eq. Q using the present experimen- 
tal velocity field measurements. The drag coefficient 
Cd = 2FjpU^D shown on Fig. [s] (top) was determined 
for the volume around the cylinder limited hy x = —2D, 
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FIG. 5: (Color) Top: Drag coefficient for // = 0.69 for 
Bernoulli-NS evaluation of pressure and for the 'flux' Equa- 
tion Q. Cdo corresponds to the non- forced flow and fit curves 
show the general behavior for this case. Bottom: Drag coef- 
ficient for // = 0.69 and for // = 2.40. 
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FIG. 6: (Color) Spectrograms for (Top) // = 0.69/o and (Bot- 
tom) = 2, 40/o. The power density spectra are represented 
as colour levels (in a log scale) on a frequency content, forcing 
amplitude (nonlinear scale) map. 



X = and y = ±41). For the non- forced case, the ref- 
erence value Cdq = 1-53 is in reasonable agreement with 
the literature (confinement effects on the Uq value have 
been taken into account in the calculation of Cdq)- The 
error was estimated from the rms value of for a varia- 
tion of ID on the control volume boundaries. The error is 
slightly lower for the impulse equation, so we adopted it 
to estimate the drag on the cylinder for different forcing 
parameters. 

Fig. [5] (bottom) shows the evolution of the drag force 
as we increase the forcing amplitude A for two cases 
// = 0.69/o and // = 2.40/o, representative of forcing 
below and above the natural frequency, respectively. For 
frequencies // < /o the drag force increases strongly af- 
ter the lock-in threshold. Qualitatively different, for />/o 
the force decreases only when we are sufficiently near the 
threshold, attaining its minimum not far passed from Ac- 
The maximum drag reduction of around 20% is consis- 
tent with what has been reported in other experimental 
works p!9l I23] as well as numerical simulations fT3| [26 j . 
Comparing the drag chart in the {ff^A) space obtained 
by [26j from numerical simulations with Fig. [2] suggests 
that for ff > /o the drag minimizes for regions close to 
the lock-in threshold, where the global fluctuations are 



reduced. On the other hand, maximum drag is obtained 
well inside the lock-in region for ff < /o as in [26j (see 
Fig.[2|. 



V. SPECTRAL ANALYSIS OF THE FORCED 
WAKE 

The lock-in threshold can be explored by inspecting 
the power spectral density on the flow domain. For each 
point in the parameter space we calculate the spectrum 
of Ux for every point in the domain. The global spectrum 
is then obtained as the sum of all points spectra. In addi- 
tion to this global spectrum that contains all frequencies 
present in the flow, we define two "spectral probes" in 
two regions chosen considering the spatial development 
of the flow: one in the cylinder vicinity, and the second 
one in the location of the maximum strength of velocity 
fluctuations. A synthesis of all the global spectra is repre- 
sented in the spectrograms shown in Fig. [6] for the same 
two forcing frequencies that have been analyzed in the 
previous sections. The spatial modification of the wake 
structure caused by the forcing is presented alongside the 
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(a) A = 0.10 




(b) A = 0.40 




(c) A = 3.00 

FIG. 7: (Color) Left: Power spectrum for a forcing frequency 
of 0.69/o for three forcing amplitudes. Black Line: Global 
spectrum; Blue line: spectrum around the global mode maxi- 
mum; Red: Spectrum in the cylinder vicinity. Right: Contour 
levels of rms{ux) with the positions of the spectral probes 
marked in blue and red. The parameters of the forcing for 
each case are indicated in the lock-in diagram. 
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(c) A = 4.00 

FIG. 8: (Color) Left: Power spectrum for a forcing frequency 
of 2.40/o for three forcing amplitudes. Black Line: Global 
spectrum; Blue line: spectrum around the global mode maxi- 
mum; Red: Spectrum in the cylinder vicinity. Right: Contour 
levels of rms{ux) with the positions of the spectral probes 
marked in blue and red. The parameters of the forcing for 
each case are indicated in the lock-in diagram. 



corresponding spectra in Figs. [7| and [Sj For each forcing 
amplitude the root mean square of Ux is displayed by its 
contour levels. 

Fig. |6ja) corresponds to a forcing frequency of 0.69/o 
and represents the characteristic behaviour for ff < /q. 
The first rank {A = 0) is the non-forced case where only 
the peaks on the natural frequency and its harmonics are 
observed. Next, under forcing, even for amplitudes a low 
as A = 0.05, the spectrum is modified by the appearance 
of the forcing frequency, its harmonics and linear com- 
binations of forcing and natural frequencies, kifj + /C2/0 
(where ki and /c2 are integers). Fig. 7 (a) [ shows a typical 
power spectral density for this regime. We observe also in 
Fig. [6]^a) that the natural frequency value increases when 
the forcing amplitude approaches the lock-in threshold. 
The peak for /o decreases its intensity as the forcing am- 
plitude increases, until the forcing frequency peak and 
its harmonics are the only present on the flow {A > 0.3 
for this case). As the global mode vanishes, the flow is 
locked on the forcing frequency (see Fig. |7(b)'] ). The peak 
on the forcing frequency continues to grow and its am- 



plitude reaches about 10 times the non forced case. For 
higher amplitudes {A > 2 f or thi s case) the spectrum be- 
comes continuous (see Fig. |7(c) ), with the forcing peaks 
distributed over a slope that clearly links higher energetic 
scales (low frequencies) to smaller scales (high frequen- 
cies). A —2 slope fits the curve, as it occurs in flows 
under strong rotational forcing [35j, unlike isotropic tur- 
bulence that fits with a —5/3 slope. This state is attained 
at high amplitudes for all forcing frequencies // < /o for 
high amplitudes and is, to our knowledge, a novel ob- 
servation in the context of open flows. The spectrogram 
in Fig. [oj^b) represents the spectrum variation for a forc- 
ing frequency // = 2.40/o. Again, for amplitudes lower 
than those near the lock-in threshold, the flow presents a 
spectrum with discrete peaks kiff-\-f2fo (Fig. 8(a) ). The 
natural frequency peak decreases its energy as the forc- 
ing amplitude increases. When approaching the thresh- 
old {A > 1.1 for this case) the natural frequency di- 
minishes its value to O.88/0. This has been discussed 
by [Ij as a consequence of the wake re-stabilization that 
changes the linear global frequency selected by the base 
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flow. Above the lock- in threshold, like on Fig. |8(b)| for 
A = 1.60 > Ac, the flow contains only the forcing fre- 
quency and its harmonics, and energy level of the fluc- 
tuations is much lower, about 10 times smaller than the 
corresponding for the non forced case. Another thresh- 
old appears as f or A > 3.2 where the spectrum becomes 
continuous (Fig. 8(c)| ), but this is owing to strong fluctu- 
ations that come from the far wake, unlike the precedent 
case where they were generated directly by the forcing 
so the phenomenon is qualitatively different. We notice 
from the spectrogram of Fig. |8(b)] that a peak on ~ ///2 
flrst emerges prior to the transition to a continuous spec- 
trum for higher amplitudes. 

It is worth mentioning that similar dynamics have been 
reported in the case of enclosed swirling flows under har- 
monic forcing [36j as well as on other forced systems [37J . 
Indeed, the quasiperiodic behavior for ff and /o that 
characterizes the transition to lock-in has been reported 
by [36j , who have described in detail the period doubling 
process. Further investigations on Re < Rcc may shed 
more light on the mechanisms of the transition process 
and the similarities of enclosed swirling flow studied by 



VI. CONCLUSIONS 

Spectral analysis was shown to be a useful tool for the 
analysis of PIV experimental data from a forced cylinder 
wake. Lower Reynolds number than previous works, al- 



low us to characterize the flow's critical behavior more 
accurately. In addition to reflning previous studies on 
the global mode properties of forced wakes and revisit- 
ing their scaling properties, the present results have in 
particular allowed to shed new light on previously un- 
explored phenomena related to the transitions between 
globally unstable and locked-in states in the parameter 
space of the forcing (/, A). The appearance of a contin- 
uous spectrum for large forcing amplitudes was observed 
in two different situations with presumably two different 
physical explanations: for / < /o, the large forcing am- 
plitude determines each vortex shed in the near wake to 
be very intense and thus to be destabilized, split and mix 
rapidly giving rise to a turbulence-like window in the pa- 
rameter space. This regime is very easy to obtain and 
can be used extensively to generate turbulent behavior 
at moderate Reynolds numbers. On the other hand, for 
/ > /o, passed the lock-in threshold, under increasing 
forcing amplitude, the flow is subject to large fluctua- 
tions coming from the far wake. 

Finally, we have evaluated the drag force from the ve- 
locity fleld bringing experimental evidence to suggest a 
relationship between the drag minimum and the lock-in 
threshold. 
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